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PARTIAL DATA INVERSE PROBLEMS FOR MAXWELL EQUATIONS 

VIA CARLEMAN ESTIMATES 

FRANCIS J. CHUNG, PETRI OLA, MIKKO SALO, AND LEO TZOU 


Abstract. In this article we consider an inverse boundary value problem for the time- 
harmonic Maxwell equations. We show that the electromagnetic material parameters are 
determined by boundary measurements where part of the boundary data is measured 
on a possibly very small set. This is an extension of earlier scalar results of Bukhgeim- 
Uhlmann and Kenig-Sjostrand-Uhlmann to the Maxwell system. The main contribution 
is to show that the Carleman estimate approach to scalar partial data inverse problems 
introduced in those works can be carried over to the Maxwell system. 


1. Introduction 


In this paper we discuss an inverse problem for the time-harmonic Maxwell equations 
with partial data, and show that the electromagnetic material parameters are determined 
by measurements on certain parts of the boundary. The result is new even for the case 
of bounded domains in but it will be convenient to formulate it more generally on 
compact manifolds with boundary. 

Let (M, g) be a compact oriented Riemannian 3-manifold with C°° boundary. The 
electric and magnetic helds on M are described, respectively, by 1-forms E and H which 
satisfy the Maxwell equations in M: 


( 1 . 1 ) 


J *dE = iujjdH, 

I *dH = —iooeE. 


Here a; > 0 is a hxed frequency, d is the exterior derivative, and * is the Hodge star 
operator for the metric g. The material parameters £ and p are assumed to be complex 
valued functions in C^[M) and to satisfy 

Re(e) > 0, Re(/i) > 0. 


The inverse problem is formulated in terms of partial measurements of the boundary 
tangential traces tE and tEl of E and H. Here, the tangential trace of a fc-form rj is 
dehned by 

t : r] i*ri, 

where i : dM —)■ M is the inclusion map. If Ti and r 2 are open subsets of dM, we dehne 
the partial Cauchy data set 

= {(fE|r^, tRlpa); {E,H) G H‘^{M,A^My solves fll.ll) and supp(fE) C Ti}. 
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As described in the end of Section [2l if a; is outside a discrete set of resonant frequencies, 
then the knowledge of is equivalent with knowing the partial admittance map 

^ tif |r„ 


where / is a 1-form on DM with supp(/) C hi, and (i?, H) is the unique solution of fll.ll) 
satisfying tE = f on dM. This corresponds to prescribing tE on Ti and measuring tH 
on r 2 . The inverse problem is to determine the coefficients e and p from the knowledge 
of for some choices of Tj. 

One may think of the above inverse problem for Maxwell equations as a generalization 
to systems of the inverse conductivity problem introduced by Calderon [CaSOj . In the 
case where M is a bounded domain in and g is the Euclidean metric, it was proved 
in [OPS93j that measurements on the full boundary (Ti = r 2 = dM) determine £ and 
fi uniquely. Earlier results include |SIC92t ISU92t ICP92j , and a simplihed proof was pre¬ 
sented in [OS96j . Stability results for this inverse problem are in |CalOt ICallj . boundary 
determination results are in [Mc97[[JM00j . the case of chiral media is considered in |Mc00j . 
and a recent uniqueness result for coefficients is given in |CZ13j . 

The inverse problem for Maxwell equations on manifolds was discussed in [QPS03] . The 
hrst uniqueness results in the non-Euclidean case were given in [KSU09j . again for the 
full data case Pi = r 2 = dM, when M is an admissible manifold in the following sense. 


Definition 1.1. A compact Riemannian manifold {M,g) with smooth boundary is called 
admissible if (M, g) is embedded in (T, g) where T = M x Mq and g = c{e®go), where c is 
a smooth positive function, (M, e) is the Euclidean line, and (Mo,5'o) is a simple manifold 
(a compact manifold with smooth strictly convex boundary such that the exponential map 
expp is a diffeomorphism onto Mq for each p G Mq). 

Locally, admissibility implies that in some local coordinates x = {xi,x') the metric 
takes the form 

where gQ{x') is some positive dehnite (n — 1) x (n — 1) matrix. Conversely, any metric 
that takes this form in some local coordinates is admissible provided that gQ is simple. 
Admissible manifolds include compact submanifolds of Euclidean space, hyperbolic space, 
and minus a point. They also include sufficiently small submanifolds of conformally 
flat manifolds and warped products. For more information see [DKSU091 lAFGRldj . 

In this article we will improve the results of |OPS93[ IKSU09j by considering the partial 
data problem where Pi = dM but r 2 is an appropriately chosen open subset of dM which 
is dehned by a limiting Carleman weight (LCW) of the manifold M. The notion of LCWs 
was introduced in connection with the Calderon problem with partial data in |KSU07j 
in Euclidean space. LCWs were analysed in detail in [DKSU09] also on manifolds. In 
particular, the class of admissible manifolds emerged in |DKSU09] as a natural class 
where LCWs exist and one may expect to be able to solve related inverse problems. Note 
that any admissible manifold has global coordinates {xi,x') where xi is the Euclidean 
direction. It was shown in |DKSU09] that (p{x) = xi is a natural LCW in this setting, 
and in this paper we will always assume that the LCW is given by (p{x) = Xi. 

We dehne the “front face” of the boundary with respect to the LCW (p by 

■= {P e dM I {dp{p), v) > 0}. 
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Here v is the 1-form corresponding to the unit outer normal of dfl. Our main result is 
the following theorem. 


Theorem 1.1. Let {M,g) be an admissible 3-manifold and (p{x) = Xi a LOW on M. Let 
Ej, fij G C^{M) be complex functions with positive real part, and let a; > 0. Suppose that 
the corresponding Cauchy data sets satisfy 

^dM,F _ f^dM,F 
Ol — 0-2 

for some open neighborhood F of in dM. Assume in addition that Ei = E 2 and fii = /i 2 
to second order on dM. Then Ei = E 2 and pi = /12 in M. 

We mention a particular case of the above theorem in (ch(f2) denotes convex hull). 

Theorem 1.2. Let C be a bounded domain with (7°° boundary, let Ej,g.j G 
be complex functions with positive real part, and let ca > 0. Given xq G \ch(f2), dehne 

F{xo) = {x G ; (x — xo) ■ p(x) < 0}. 

Suppose that the corresponding Cauchy data sets satisfy 

^dn,F _ ^dn,F 
L-i — O2 

for some open neighborhood F of F{xo) in dfl. Assume in addition that Ei = E2 and 
hi = h 2 to second order on dfl. Then Ei = E 2 and pi = /i 2 in fl. 


This result involves the logarithmic Carleman weight in and corresponds to one of 
the results of [KSUOTj for the partial data problem for the Schrodinger equation. Note 
that if is strictly convex, then any open subset of dfl can serve as F above. 

A standard method to study such partial data problems is to apply the idea developed 
in [BU021IKSU07] for scalar equations, where Carleman estimates with suitable positive 
boundary terms were used to suppress the unknown information. The partial data inverse 
problem for Maxwell equations presents several challenges when one tries to apply these 
methods directly. The hrst difference is the fact that the principal part of this equation is 
a hrst order system rather than a second order scalar equation, and Carleman estimates 
with boundary terms for systems seem to be more involved than those for scalar equations 
(see [E1081 ISTOQj for some such estimates). Secondly, the fact that ellipticity is somewhat 
’’hidden” in the system makes it difficult to obtain a estimate similar to the one used in 
|KSU07] which could be directly applied to the Maxwell system. 

We circumvent these difficulties by hrst using ideas from [OS961IQPS031 IKSUOQj to show 
that the problem reduces to constructing suitable solutions to a Hodge Dirac system 


(P -|- W)u — 0, ~ 0- 

We will then apply the recent work |CST13j where the authors studied Carleman estimates 
and complex geometrical optics solutions for the Hodge Laplace operator with relative 
and absolute boundary conditions. The crucial point is that the Carleman estimate with 
boundary terms proved in [CST13] for the Hodge Laplacian is sufficiently powerful to yield 
information also in the Maxwell case. However, we stress that the present situation for 
Maxwell equations does not reduce to that of |CST13j , since the boundary measurements 
for Maxwell determine in a sense only part of the boundary measurements for the Hodge 
system. It is essential to use the special structure of the Maxwell system to show that 
the coefficients are uniquely determined. 
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The main contribution of this work is to show that the Carleman estimate approach 
to scalar partial data inverse problems introduced in |BU021 IKSUOTj can be carried over 
to the Maxwell system. A Carleman estimate approach for a different first order system, 
related to the Pauli Dirac operator, was presented in |ST09j where the method involved 
decoupling the Pauli Dirac operator into a second order differential operator with the 
Laplacian as its principal part. A suitable Carleman estimate was then applied to the 
decoupled equation to recover the coefficients. The partial data problem for Maxwell 
equations was also studied in |COS09j in the case when M is Euclidean and the inaccessible 
part of the boundary is a portion of a hyperplane or a sphere, following the scalar approach 
of pi7] . Due to the partial symmetry of the domain one can reflect across the flat part 
of the boundary and reduce the problem to a full data problem. Partial data results 
are also known for a two-dimensional Maxwell system [lY 14aj and for Maxwell equations 
in a waveguide |IY 14b] . based on the two-dimensional partial data results of |IUY10] . 
and in the case where the parameters are known near the boundary jBMR14j . Recently, 
extensions of the Carleman estimate and reflection approaches for partial data problems 
were introduced in [KS13j and |IY13j : some of these extensions for the Maxwell system 
were considered in |IY14cj . We also remark that the methods in the recent paper |DKLS13] 
might allow to relax to some extent the admissibility assumption in Theorem 11.11 
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Finland (Centre of Excellence in Inverse Problems Research), and F.C. and M.S. were 
supported by an ERC Starting Grant (grant agreement no 307023). L.T. was partly sup¬ 
ported by the Academy of Finland (decision no 271929), Vetenskapsradet (decision no 
2012-3782), and Australian Research Council Future Fellowship (fellowship no G161397). 
F.C. would like to acknowledge the University of Jyvaskyla for its hospitality on subse¬ 
quent visits. 


2. Reduction to the Hodge Laplacian 


In this section we employ ideas from |QS961 IOPS031IKSU09] to show that the problem 
reduces to constructing suitable solutions of the Hodge Dirac and Schrodinger operator. 
The notation follows closely |KSU09j and we refer to that article for more details. If E and 
H are complex 1-forms in M, we consider the graded form X = ^ + E + *H + *4/ where 
$ and 4/ are complex scalar functions. We write graded forms in the vector notation 

Y = ( $ I *4/ E ) . 


The line denotes that the forms of even order and odd order are grouped together, which 
will result in a block structure for the equation. We define the following Hodge Dirac 
operator acting on graded forms in M, 


P=l(d-6) 

i 



( 

-6 \ 

1 


-6 d 

i 

d 



y d —6 

/ 


The Dirac equation which is most closely related to Maxwell is 


{P + V)X = 0 
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where V is the operator acting on graded forms as the matrix operator 


/ —OJfi 

*Da A * 

—UJfl 

*Da A * 

P/3A 

—uje 

\ D(5^ 

—bJS 


and where D = jd, a = log e, and fd = log /r. One easily sees that 


r (P + V)X = 0, r *dE = iufiH, 

\ X = {0 *H \*0 E y ^ \ *dH = -iujsE. 


The advantage here is that P is an elliptic operator. 

For the reduction to the Schrodinger equation with Hodge Laplacian, we consider 
rescaled solutions 


Y 



and henceforth relate X and Y in this way. We will write 


r = ( yo y2 I y3 yl ^ 

where is the fc-form part of the graded form Y. It is easy to see that 


(P + 1/)A = 0 ^ 

(P + 1+)F = 0 

where W is the potential 




/ 

*Da A * \ 

1 

+ 

to 1 


*Da A * —DaA 

D/dA 



y D/dA *D(d A * 

/ 


and A = 

The reason for using the rescaled solutions Y is that there is a good reduction of the 
Dirac equation {P + W)Y = 0 to a Schrodinger equation with no hrst order terms (see 
jKSUnQl Lemma 3.1]). 


Lemma 2.1. Denote by —A = dS + Sd the Hodge Laplacian acting on graded forms. We 
have 


{P + W){P -W^) = -A + Q, 
(P - W%P + W) = -A + Q', 
{P + W*){P-W) = -A + Q, 














6 


CHUNG, OLA, SALO, AND TZOU 


where Q, Q', and Q are continuous potentials (endomorphisms of AM) with 



/ Aa + \ {da, da) 

0 

• 


2 1 
"A 

0 

• 

• 

• 

• 

• 

A/3 -f i(d/3, d/3) 

0 


\ 

• 

0 

•/ 


/ A/3 — ^(d/3,d/3) 

0 

0 

0 \ 

2 

• 

• 

• 

• 

0 

0 

Aa — i(da, da) 

0 


\ 

• 

• 

•/ 


and • denote hounded coefficients involving derivatives of e and /i up to second order. 


Remarks. The coefficient functions denoted simply by a “bullet” above are not needed 
in our argument. They are easy to compute though if one so wishes. Also observe that 
due to this lemma, if Z is a solution of {—/A + Q)Z = 0, then Y := {P — W)Z would solve 
{P + W*)Y = 0. Finally, the precise form of Q is not needed in this work, so we will not 
bother the reader by writing it out explicitly. 


The above arguments can be used to reduce the Maxwell equations to the Hodge Dirac 
and Schrodinger equations. We now want to give a similar reduction on the level of 
boundary measurements. In fact, we will prove an integral identity showing - see Lemma 
12.31 - that if the partial Cauchy data for two sets of coefficients coincide, then the dif¬ 
ference of the corresponding potentials Qj is orthogonal to a product of solutions to the 
Schrodinger and Dirac equations. This result is a version of |KSU091 Lemma 3.2] adapted 
to the partial data problem. 

If U and V are graded forms on M, we use the notations 

„ 3 r ^ 

UA)= / E<'^tU)£iV', (U\V)aM= 'T(U\V’)dS. 

JM JdM 

for the L^-inner products. We also use the same notation in the case when U and V 
are in appropriate dual Sobolev spaces. Note that in the second expression, U and V do 
not have be graded forms of the manifold dM (for instance, the 1-form parts may have a 
normal component). We have the integration by parts formulas f |Ta991 Section 2.10]) 

(2.2) {dU\V) = {uA U\V)9M + {U\6V), 

(2.3) {6U\V) = -{iuU\V)9M + {U\dV) 

where u is the outward normal and i^ is the contraction with u. Also, we define 

HA = {ue L2(M,AM); Au G L‘^{M,AM)} 


and equip this with the natural graph norm 

II«IIL = 11“ 


+ ||A«|p. 


Here and in the rest of this paper, notation of the form || ■ || and || ■ \\m should be taken 
to indicate the relevant norm unless otherwise specihed. 

One feature of the method in this paper is that we need to deal with pairs of solutions, 
where one solution can be very nonsmooth (only in general). For this reason we will 
need the following result concerning smooth approximation and traces in the Ha space. 
The proof is given in the appendix. 
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Lemma 2 . 2 . The set AM) is dense in H/^. The trace maps u i—)■ {tu,t6u), 

M I—)■ (t*u,t6 *u), M I—)■ (tu, tdu) and m ha- (t * u,td* u), initially defined on AM), 

have unique hounded extensions as maps 

Ha H-^/‘^{dM,AdM) x H-^/^{dM,AdM). 

Given a pair of electromagnetic parameters {Sj, pj}, j = 1, 2, we define the correspond¬ 
ing matrix potentials Vj, Qj and Wj be replacing e: and fv respectively with Sj and fVj in 
dehnitions above. The next result is the main integral identity. 

Lemma 2 . 3 . Let G C^{M), j = 1, 2 , be complex functions with positive real parts, 
and assume that 

(2.4) Si = 82, /ii = /i2 to second order on dM. 

Let T2 be any open subset of dM, and assume that . Then 

{{Qi-Q2 )Zi\Y2 ) = 0 

for any graded forms Zi G and I2 £ such that 

(—A -|- Qi)Zi = 0, 

Yi = {P- Wl)Zi satisfies Yf = Yf = 0, 

and 

Y2 = {P — 1 T 2)^2 for some Z2 G Ha with (—A -|- < 52)^2 = 0, 

supp(ty2) c r2. 

Remark. By the assumptions for Sj and p,j, we have that Qi — Q 2 is and vanishes 
on dM. Thus {Qi — Q2)Zi is H^ and vanishes on dM, and the pairing ((Qi — Q2)Zi\Y2) 
where Y 2 G makes sense. Note also that since Z 2 G Ha, Lemma [2.21 guarantees that 
the boundary value tY2 exists as an element of H~^/‘^. Here, Zi and Yi need to be related 
to a solution for Maxwell, but Y2 only needs to solve a Dirac equation. This simplihes 
the recovery of coefficients. 

Proof of Lemma By the discussion earlier in this section, if Yi is as stated, then 
Yi G H^ solves (P + lTi)yi = 0 and gives rise to a H^ solution Xi = ( 0 *Hi | 0 Pi 
of the Maxwell equations (P -|- Vi)Xi = 0. By the assumption on Cauchy data sets, there 
is X 2 = ( 0 *H 2 I 0 E 2 y & H"^ satisfying 

{P + V 2 )X 2 = 0 , 

tE2\dM = tEi\QM, tH2\r2 = ^hfi|r2- 

Let Y 2 G H"^ be the solution of (P -|- 1 ^ 2)^2 = 0 corresponding to X 2 . 

We hrst claim that the assumption on Cauchy data sets implies 

(2.5) ((1Hi-1H2)Fi|F2)m = 0. 

Here, 1^2 is only in H~^ since Y2 = {P — W2)Z2 where Z2 G Ha, but the pairing is well 
dehned because the quantity {Wi — ^ 2 ) 1^1 is in H^ with vanishing boundary values since 
El = 82 and Pi = p 2 to second order on dM. To deal with the nonsmooth solution Y 2 , 
we use Lemma [2.21 and choose Zi G C°°{M,AM) with Zi —)■ Z 2 in Ha as I —>■ 00 , which 
implies in particular that (P — 1 ^ 2 )-^/ Y 2 in H~^. Noting that 

(P + W2){Yi - Y2) = (W2 - Wi)Fi, 
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we compute 

(( 1^2 - Wi)Yi\Y2)m = \im{{W2 - Wi)Yi\{P - W2 )Zi)m 

= lim((P + W 2 ){Y^ - Y 2 )\{P - W 2 )Zi)m 


= lim 

i 


-((z/ A +*,)(Fi - Y2)\{P - W2)Zi)sm + (n - mP + W;){P - W2)Zi) 


M 


In the last line we used fl2.2p and fl2.3p . For the last term, we have 

( 2 . 6 ) (P + iy2*)(P - 1 ^ 2 )^; = (-A + Q2)Zi ^ (-A + <52)^2 = 0 in 

so the second term in brackets goes to zero as / —)■ cxo. For the boundary term, we claim 
that 


(2.7) t{Yi-Y 2 ) = Oon^M 

(2.8) t * (Fi — 1 ^ 2 ) = 0 on F 2 

(2.9) tY 2 = 0 on F 2 . 

Assuming these, one has u A (Yi — Y 2 )\dM = 0 and P(Fi — F 2 )|r 2 = 0, and the boundary 
term becomes 

\im-{Y,-Y2\iy A{P -W2)Zi)r:^. 

i i ^ 

Here t(P — W 2 )Zi\y<^ -a tF 2 |r| = 0 in py Lemma [2.21 Since Fi — Y 2 \dM is in 

the limit of the pairing on F 2 is well dehned and in fact zero. This proves (12.51) modulo 
(I2.7p - (l2.9p . To check (12.7p . we see that 

z/A (Fi-F2) = ( 0 e\'\^{El-E2)\*^^T{^.Hl-H2) 0 )• 

From the hypotheses, z/ A (Pi — p 2 )|aM = 0, and using the surface divergence as in the 
proof of [KSU09, Lemma 3.2], we also have (z^. Pi — P 2 )| 9 Ar = 0. Therefore (12.7p is 
satished. Meanwhile, 

z/A*(Fi-F2) = ( 0 /xJ/Va (Pi-P2) I *£}/^(z/,Pi-P2) 0 ). 

From the hypotheses, z/ A (Pi — P 2 )|r 2 = 0, and using the surface divergence as in the 
proof of [KSU09, Lemma 3.2], we also have (z/. Pi —P 2 ) |r 2 = 0. Therefore (12.81) is satished. 
The equation (12.91) follows straight from the dehnition of F 2 . 

Having proved (12.5p . we next claim that 

(2.10) ((IFi - W2)Y^\Y2)m = ((Qi - Q2)ZAY2)m. 

Note that the pairing on the right is well dehned, since {Qi — Q 2 )Zi is in P^ with vanishing 
boundary values. Since Qj = Wj o P — P o Wj — WjWj, we have 

(IFi-IF2)Fi = (1Fi-IF2)(P-IF*)Pi 

= (Qi + P o IF()Pi -{Q 2 + P 0 + W2{W^ - IF())Pi 

= (Qi - Q2 )Zi + (P + W2)WlZi - (P + H^2)kF2%. 
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Choosing again Zi G AM) with Zi —)• Z 2 in Ma, we have 

- W2)Yi\Y2)m = - W2)Yi\{P - W2)Zi)m 

= lim [((Qi - Q 2 )Z,\{P - W 2 )Zi)m + {{P + W 2 ){Wl - W^)Z,\{P - W 2 )Zi)m] • 

The hrst term in brackets has limit ((Qi — Q 2 )Zi\Y 2 )m- After integrating by parts using 
(I2.2p and (I2.3p . the second term becomes 

-((p A - W^)Zi\{P - W2)Zi)aM + ml - Wl)Z,\m + Q2 )Zi)m • 

i 

Since at dM to hrst order, the boundary term vanishes. Also, by fl2.6p the 

second term has limit zero. This proves fl2.10p . and combining this with (12.51) proves the 
lemma. □ 

We now explain the connection between the Cauchy data set and the admittance map, 
which is a standard way for expressing boundary measurements for Maxwell equations. 
First we introduce some function spaces. The surface divergence of / G H^{dM, A^dM) 
is given by 

Div(/) = {ddMf, dS) 
where dS is the volume form on dM. Dehne the spaces 

Him) = e H\M,A^M) ; Dw{tu) G 

THlilidM) = {/ G H^/\dM,A^dM ); Div(/) G H^/\dM)}. 

By [KSUOQl Theorem A.l and subsequent remark], there is a discrete set S C C such that 
if (V G C \ S, then for any / G THl{l{dM) the Maxwell equations (11.11) have a unique 
solution (E, H) G Hl^^{M) x Hl^^{M) satisfying tE = f on dM. Assume that oj is not a 
resonant frequency, and let Ti and r 2 be nonempty open subsets of dM. If THI{1^{Ti) 

consists of those elements in THl{l{dM) that are compactly supported in Ti, we may 
dehne the partial admittance map 

Ar,.r, . ^ tH\r, 

where [E, H) is the unique solution of (11.11) satisfying tE = f on dM. Note also that any 
solution (E, H) G X H‘^ of (11.11) is in x by using the surface divergence as 

in |KSU09l Lemma 3.2]. 

Now if bj is not a resonant frequency, it is clear that is the graph of 

= {(/, ; / e rffD(lc(ri)}. 

Thus and A^^^^ are equivalent information in this case. 

3. Carleman estimates and CGO solutions 

Let (M, g) be an admissible 3-manifold as in Dehnition ll.il so that M CC M x Mq and 
g = c(e © go) where {Mo,go) is simple. Let {xi,x') be global coordinates in R x Mq, and 
let ip{xi,x') = xi be the natural LCW on M. 

It is convenient to begin with a reduction to the case where c = 1. This can be done 
using the following lemma from |KSU09j . 


lim 

i 























10 


CHUNG, OLA, SALO, AND TZOU 


Lemma 3.1. (^ [KSUOQj . Lemma 7.1) Let {M,g) be a eompact Riemannian 3-manifold 
with smooth boundary, and let e and /i be functions on M with positive real part. Let 
c be any positive smooth function on M, and let represent the Cauchy data for e, fv 

with respect to the metric g. Then 

/^Ti,r2 _^ri,r2 

Proof. Lemma 7.1 in |KSU09j is actually phrased in terms of the admittance map, which 
requires the additional condition that solutions to the Maxwell equations are unique. 
However, the same proof works when the lemma is phrased in terms of the Cauchy data. 
The key fact is that for a fc-form u, 

*cgU = U, 

SO {E,H) satishes fll.ip with metric eg and coefficients £ and /i if and only if {E,H) 
satishes fll.ip with metric g and coefficients and □ 

Thus it suffices to prove Theorem 11.11 in the case where c = 1. We will from now on 
assume that g = e® g^ where {Mo, go) is a simple 2-manifold, and that we are working 
with the LCW (p{x) = Xi. We will also assume that (Mo,5'o) is another simple 2-manifold 
which is slightly larger than Mq. 

As indicated in Lemma [2.31 we will need to construct two types of complex geometrical 
optics solutions. First we need sufficiently regular solutions to the Maxwell equations, with 
no restrictions on their boundary values. By Section [2] it suffices to construct solutions Z 
to the equation (—A -|- Q)Z = 0 in M such that Y = {P — W*)Z satishes = 0. 

Such solutions were obtained in [KSUOQl Theorem 6.1a]. Second, we need solutions of the 
Dirac equation (P -|- W*)Y = 0 satisfying tY\r = 0. Here T C dM is meant to be a open 
neighbourhood of dM \ F, where F is as in the statement of Theorem 11.11 By proper 
choice of T, we can ensure that 

(3.1) {dip, I/) < 0 on T. 

The next theorem states the existence and basic properties of both types of complex 
geometrical optics solutions. 


Theorem 3 . 1 . Let £,/i G C^{M) have positive real parts, and let W and Q be as in 
Section m Fix p G Mq \ Mq, and let {r,9) be polar normal coordinates in {Mo,go) with 
center p, so that {xi,r,9) are global coordinates near M. Let so,to he real constants, let 
A > 0, and let h{9) G C'°°(5i). 

(a) For |r| sufficiently large and outside a countable subset of M, one can construct 
solutions to the equation (—A -|- Q)Z = 0 in M satisfying 

(P + W)F = 0, Y = {P -W^)Z, 

yo = ^ 0 


where Z G H^{M, AM) has the form 
(3.2) Z = 9)\-^/^P^^^^+"'>b{9) 

with IIPIIhs < for 0 < s < 2 . 


f So \ 
0 


to * 1 

V 0 / 
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(3.3) 


(b) Let r C dM be an open set so that (13.11) holds for ^{x) = xi. For r > 0 sufficiently 
large, there exists a solution Z G i^A to (—A + Q)Z = 0 in M so that 

{P + W*)Y = Q, Y = {P-W)Z, 

and Y has the form 

/ So \ 

-isodx^ A dr 




fn * 1 


+ e 




{R + R') 


\ito * dx^ A drJ 

with ||.R||j:/-i < Ct~^R, ||i?'||L 2 < Ct~^R, and Y satishes the boundary condition 

fF|r = 0. 


Proof of Theorem \3.1\f a). It is proved in |KSU09( Theorem 6.1a] that there is a solution 
Z G AM) of this form. That theorem only states bounds for R, but the proof 

of the theorem relies on |KSU091 Proposition 5.1] which actually gives the if* bounds for 
0 < s < 2. It remains to show that under the assumption e,/i G C^{M) (which implies 
Q E C^) one has in fact Z G H^{M, AM). This follows from interior elliptic regularity 
since the solution in |KSU091 Theorem 6.1a] is constructed in a neighbourhood of M, and 
Z satishes —AZ = —QZ where QZ G in this neighbourhood. □ 


In the remainder of this section we will prove Theorem I3.1f b). This will be done by 
constructing a solution for the equation (—A + Q)Z = 0 of the form 

(3.4) Z = e^^^{A + R), tZ\r = t{S + iW)Z\r = 0, 

for some appropriate amplitude A chosen so that 

e-™H-A + Q)e™M = Oi2(M)(l). 

Then Y = {P — W)Z solves {P + W*)Y = 0 and satishes the desired boundary conditions. 

The remainder term R in 03.41) will be constructed using the following Carleman esti¬ 
mate for graded forms with suitable boundary conditions. If m is a graded form we write 
u± = u A Ru and u\\ = u — u±. The proof of the following estimate was hrst given in 
|CST13j . Here || • || = || • ||m and || • \\qm are the relevant norms in M and dM. 

Proposition 3.2. ( [CST131 Theorem 7.2]) Let {M,g) be an admissible manifold as in 
Definition \TT\ and (f{x) = Xi a LOW. Let Q and a be L°° endomorphisms on the space 
of graded forms. There is C > 0 such that for all smooth graded forms u satisfying 

-|- tau±\r = fwlr = 0, 


u\dM\r — '^uu\oM\r — 0, 

and for all r > C, we have the following estimate 

r||M|| -L ||Vu|| +T^^‘^\\ti^u\\dM + V^llV'tiuulldM + V^lltVuUiilldM < C||e^‘^(-A-hQ)e“^‘^M||. 

Proposition 13.21 relies on the fact that (13.Ih holds on P - this is one of the key reasons 
for choosing F as we did in Theorem 11.11 




























12 


CHUNG, OLA, SALO, AND TZOU 


3.1. Construction of the amplitude. If Q is an endomorphism of AM, we look 
for a solntion Z of 

{-A + Q)Z = 0 inM 

having the form 

Z = e^^^{A + R) 

where r > 0 is large. To this end we will construct an amplitude A such that 

Recall that since (M, g) is an admissible manifold and we have reduced to the case c = 1, 
the metric can be written globally as 5 ^ = e © 5^0 for some conformal factor c and metric 
Qq on the simple 2-dimensional manifold Mq. 

We now construct the amplitude A mentioned above. 

Proposition 3 . 3 . For any smooth functions V,bi,bl G there exist graded forms 

A = A^ + A^ + *B^ + satisfying 


||e-""iAe™M||M = 0(1), ||A||m = 0(1) 

of the form 

= \go\~^^'^e'''"''b^{e)dxi + {\go\~^^%l{e) dr + 

= \go\-^/^F^^b\e)dx^ + F^^{\go\-^/%l^e) dr + \go\^/%l{e) dO). 

Proof. Any 1-form A^ on M may be written as 

A^ =Aldx^ + {A^y 

with A\ = {A, dx^) and {A^)' is a 1-form on Mq that contains xi as a parameter. If A° is 
a 0-form and A^ is a 1-form, we have 

A,A° = {dl + A,0^°, 

A^A^ = (A3AJ) dx^ + dl{A^)' + A,/(A^)' 

where A^,/ is the Hodge Laplacian on (Mo,(7o)- 

We compute the conjugated Hodge Laplacian acting on 0-forms and 1-forms, 

= (-92 - 2r9i - - A,0A°, 

= [(-d" - 2r9i - r" - A^O^J] dx^ 

+(-V|^ - 2rVa, - - A,0(^')'- 

We now make the simplifying assumption that A° and A^ only depend on x'. Then the 
above expressions simplify to 

e-™i(-A3)e™M° = (-A,, - 

= [(-A,, - t^)AI] dx^ + {-A^, - T^){A^y. 
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It is enough to arrange that 

(3.5) - t^)A^\\m, = 0(1), ||A°||mo = 0(1), 

||(-Aa,/ - t‘^)A\\\mo = 0(1), IIAJIImo = 0(1), 

IK-A,,, - t^){A^)'\U = 0(1), ||(A^)1Uo = 0(1). 

The construction for the scalar parts A^ and Aj is identical: we look for a WKB type 
quasimode 

and note that 


(—Aa;/ — r^)A° = [r^(|(i-^|^ — l)a — ir(2Vvya + (A'0)a) — Aa] . 

In order to satisfy fl3.5p . we need -0 and a to satisfy the eikonal equation \d'ip\‘^ — l = 0, and 
the transport equation 2Vvya + {A'il!)a = 0. To hnd ^fJ, let {Ado, go) be a slightly larger 
simple manifold than {Ado, go), let p G Ado \ ^o, let {r,6) be polar normal coordinates in 
{Ado, go) with center p, and choose 


-(/^(r, 9) = r. 

In these coordinates, the metric has the form 

(^0 m{r, 6 ) ^ 

for some smooth positive function m. It follows that 

|#P = 1. 


We also have = dr and Ai/; = lclr(log|fi'o|) where \go{r, 6 )\ is the determinant of go- 
The transport equation becomes 

2dra+ ^dr{\og\go\)a = 0 


and it has the solution 

a = \go\~^^%0) 


where b can be any function in C°°{S^). Therefore, setting Ao = \go\~^^'^A'^^h^{ 6 ), the hrst 
line of (13.51) is satished. A similar choice for A\, substituting for 6°, would satisfy the 
second line of (13.bh . 

Moving on to (A^)', we look for an ansatz of the form 


{A^)' = e^^a', a' = ardr + aede. 


We have 


(—Aa,/ — t‘^){A^)' = — l)a' — ir[2Vvya^ + (A'0)a'] — Aa'] . 

Choosing the same function 'ip{r, 6 ) = r as above we have {dipl'^ = 1. Since 


Va.dr = 0, Va,.de = --^^(loglsfol) dO, 






14 


CHUNG, OLA, SALO, AND TZOU 


the transport equations for and ag become 

‘ZOfCbf (logl^oI)®r 0, 

2drag - ^dr{\og\gQ\)ag = 0. 

These have the solutions 

Or = \go\~^^%{0), ag = \go\^^%e{9) 

for any bg,br G C°°{S^). Setting A' = {\go\~^^^bl{6)dr + \go\^/%l{6)d6), the third line 

of (13.5|1 is satished. 

Now we can use the fact that * commutes with A to construct and with the 
same asymptotics, and hnish the proof. □ 

3.2. Construction of the remainder. We now construct the remainder term R with 
the appropriate boundary conditions and asymptotics in r. More precisely we need to 
construct R to satisfy 

e-™i(_A + Q)e^^^R = -e-™i(-A + Q)e™M 
with boundary conditions 

tR\r = -tA\r, R) + iWR) |r = + iWA) |r. 

The hrst step in constructing such an R is the following solvability result similar to 

|CST131 Proposition 7.3], obtained from Proposition 13.21 Here we write A^ = e’'‘^Ae“'^‘^. 

Proposition 3.4. Let Q and W he endomorphisms on M acting on graded forms. 
If T > 0 is large, for any graded forms F G Lf{M,kM) and /i ,/2 G AM) there 

exists a solution u G Lf{M,AM) n to 

(— A—ip + Q)u = F 

(3.6) tu\r = tfi\r, 

{e-^nde^^u + tWu,_)\r = + thP(/ 2 )±)|r 

and which satisfies the estimate 

(3.7) ||m||m A — ^IIT’IIm +-^lle "'‘^t5e'^‘^/2 + ^hP(/2)±||aM + \/d||^/i||9m^ • 

Note that the boundary conditions in fl3.6p make sense because u G guarantees 

that the traces in fl3.6p are well dehned; see Lemma [2.21 

Proof. Let H be the subspace of Lf{M,AM) dehned by 

H := {(—A<^ + Q*)v I V G C°°, tv = 0, e'^'^t5{e~'^'^v) = —tii,W*iyV, n|rc = VAru|r': = 0} 
Dehne on H the linear operator 


L((-A^ + Q*)v) := {v\F)m + (e"^^,de-"^n|t/i)r + h) + tW{fA±)r. 







PARTIAL DATA INVERSE PROBLEMS FOR MAXWELL 


15 


We see that L is well defined and bounded, since by Proposition 13.21 and by the fact that 
tijydv = —Vti^v for v as above |CST13i Lemma 3.4], the quantity \L{w)\ is bounded by 


c, 


7 ^ 


Vr\\tfi\\r + Vr\\iuv\ 


|e +fW(/ 2 )±||r 


r 


< -||(-A^ + Q*)n||M Ili^lU + V^WtfiWr + ^||e-"^t5(e"^/2) + thP(/2)T||r 


T 


By Hahn-Banach, there exists u E with the estimate 


|w||m < ^ (^lli^lk + V^WtfiWr + -^We-^ndie^^h) + thP(/ 2 )T||r) 


satisfying 


((-A^ + Q*)v\u) = L((-A^ + Q*)v) 

= {v\F)m + + (*.(e-"^t;)|M(e"‘"/2) + e"nW(/2)a)r 

for all V as in the definition of H. 

We consider first those v which are compactly supported in the interior of M to see 
that u solves (—A_(^ + Q)m = F. Thus in particular u G so the traces of are 

well defined by Lemma 1221 Now we can integrate by parts using fl2.2p . (I2.3p . the relation 
= —i^W*iyV, and the fact that v vanishes to first order on P'^, to see that 
on the boundary u satisfies 

{te^‘^u\iyd{e~'^‘^v))Y + {te~'^'^5{e^'^u)\ti^v)Y + {iyu\ii,W*i,jv)Y 
= ite^^ h\iud{e-^^v))r + (te-"^<5(e"^/2) + tW U 2 ) i_\ti,v)v 
which becomes 

{te^‘^u\iyd{e~'^‘^v))Y + + tWu±\tii,v)Y 

= {te^^h\i,d{e-^^v))r + f2) + tW{f2)±\U,v)r 

for all V as in the definition of H. Now for all 'ijj G with supp(ij,'0) C P, we may 

choose V G such that 

v\dM = 0 , t5{e~'^‘^v) = 0 , iyd{e~'^'^v) = iy'ip 

(See Lemma 3.3.2 of |Sc95] .l A brief computation shows that v vanishes to first order on 
P'^ and can be used in the above identity, which gives 

whenever supp(i,y'^) C P. Therefore we see that on P and 

{te-^^5{e'^^u)+tWu^\U,v)v = f2) + tW{f2)i\U,v)v 

for all V as in the definition of H. Applying an argument similar to |Sc95l Lemma 3.3.2] 
(see also Lemma l52|) again allows us to conclude that 

e-^^t5{e^^u) + tWu^ = f 2 ) + tW{f 2 )± 

on P. □ 


This gives the following existence result for complex geometrical optics solutions with 
boundary conditions. 
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Proposition 3.5. Let Q and W be endomorphisms acting on graded forms. Forr > 0 
sufficiently large there exist solutions to the eguation 

(-A + g)Z = 0 


of the form 

with boundary conditions 


Z = e^^^{A + R), 


tZ\i' — 0, t(5Z -\- iWZ |p — 0 

where A is as in Proposition 1 3 . 31 with R G L?‘r\e~'^^^H/^ satisfying the estimate ||i?|| < 


Proof. By Proposition 13.31 the norm of F := e A + Q)e'^^^A is 0(1) as r —)■ oo. 

Now apply Proposition 13.41 to obtain i? G fl solving 


e-™i(_A + Q)e^^iF 
tR\r 

t{e-^^^6e^^^R + iWR±)\r 


-F, 

-tA\r, 

-t(e-^^^6e^^^A + iWA^)\r 


with the estimates 

\\R\\m < -(\\F\\M + V^\\tA\\9M + ^\\e-^^^t6e^^^A + tiWA^ 
T V \/t 



□ 

Proof of Theorem \3 . iV b) . By Proposition 13.51 there is a solution Zq G PTa of the equation 
(—A + g)^o = 0 satisfying 

Zo = + ^^olr = t{5 + iW)ZQ\r = 0 

where the amplitude A is chosen to be 

A = R-^\g\-F%e)C 2 

with ^2 = —i{sQ dx^ + to * dx^) where Sq, to are real constants. 

We dehne Z := ^Zq, so Z will solve the same equation as Zq with the same boundary 
conditions. Then Y := {P — W)Z will solve {P + W*)Y = 0 in the sense of distributions. 
Since tZ|r = 0, we have Z\t = Z±\y and tdZ\Y = 0 by Lemma [5Tl Consequently 

tY\r = -{tdZ - t6Z - itWZ)\r = -U{5Z + iWZ^_)\r = 0. 
i i 

Finally, we observe that 

Y = -{P - W){eR^^+^^^\g\-F%{e)C + 

T 

r T T 

+ e™i-(P-W)(P). 
r 
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A direct computation (see jKSUn9( Proof of Theorem 6.1(b)]) gives the result. 

4. Recovery of coefficients 


□ 


Here we use the solutions constructed in Theorem 13.11 to recover the coefficients of the 
Maxwell equations from partial boundary measurements. The argument is essentially the 
same as in |KSU09j . 

Proof of Theorem \l.l[ Let Zi be the solution to {—A+Qi)Zi = 0 constructed in Theorem 
I3.1f al having the form 

( s, \ 

0 




to* I 

V 0 / 




Choose an open set Fi C dM such that C Fi CC F, and let T = dM \ Fi. Then 
holds, and we can take Y 2 to be the solution to (P + 1T2*)^2 = 0 constructed in Theorem 
l3.11 bL satisfying fY 2 |r = 0 and having the form 

/ So \ 

isodx^ A dr 


Y 2 = 


to * 1 


+ e"(^i-*")(P + P')- 


\—ito * dx^ A dr J 

(In fact we solve (P — Wf)Y = 0 where Y is as in Theorem l3.11 bL and we let I 2 = ^0 
Here Sq, to, X are any real constants with A > 0, b{9) is any function in C°°{S^), and r > 0 
can be taken arbitrarily large. 

By taking r 2 = P in Lemma [2.31 we have supp(fy 2 ) P r 2 . Thus we have the integral 
identity 

((Qi-Q2)^i|k^2)= [ {{Qi-Q2 )Z,\Y2)dV = 0. 

Jm 

Writing this out term by term we have 
(4.1) 



/ So \ 


/ So \ 

[ ({Qi- 

JM \ 

0 


—isodx^ A dr 

to * 1 


to*l 

\ 

'v 0 y 


\ito * dx^ A dr J 


{Q,-Q^)R\{R + R'))\gfP + 


IM 


' M 


(Qi-Q2)C^("^+''-)6(^) 

( 


( So \ 
0 


to * 1 

V 0 / 


(Qi-Q 2 )P| 


'M 


So \ 
—isodx^ A dr 


\R+ii)\af'* 


laY^o. 


to *1 

\ito * dx^ A dr J 

Here all integrals are with respect to the 3-form dV = dxidrdO. 

We will argue that all terms in (14.Ih except for the first one vanish as r —)■ cxo. Indeed, 
by Theorem 13.11 ||P||_us < CY~^, ||P||jii--i < and ||P'||l 2 < Using these 
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estimates for s = 1 and s = 0 we see that the second term is Here we use that 

El = £2 and /ii = /12 to second order on dM, to have Qi = Q 2 on dM. By the same 
estimates for s = 0 we see that the third term is and the fourth term is 0 (r“^). 

Thus taking r —)■ cxd in fl4.ip implies that 



(^{Qi — Q 2 ) 


/ So \ 
0 

to* I 


V 0 / 


/ So ^ 

—isodx^ A dr 
to* 1 

\ito * dx^ A dry 


dr dO 


0 


for all real So, to, A with A > 0, and all smooth functions b{9) G 

The rest of the proof is identical with the proof of Theorem 1.1 in |KSU09j : we give 
the argument for completeness. Let and be the elements of Qi — Q 2 , interpreted as 
a 8 X 8 matrix, which correspond to the (1, l)th and (5, 5)th elements, respectively. By 
Lemma 12.11 

Qa = ■^A(ai - ^2) + ^{dai,dai) - ^{da2, da2) - - £2/1-2), 

qp = ^A(/3i - ( 32 ) + ^{d(3i, d(3i) - ^(d/? 2 , d/ 32 ) - - £ 2 / 12 )- 


With the two choices (so,^o) = (1,0) and (so,^o) = (0,1), the special form of Qi and Q 2 
in Lemma [2.11 shows that we obtain the two identities 


Im dxi dr d9 = 0, 

Im e*^*''^^’'’*'’^&(d)g/ 3 (a;) dxi dr d9 = 0. 


Let T = M X Mo. Since Qi\dM = Q 2 \dM, the zero extensions of qa and qp to T \ M 
are continuous functions and the integrals above may be taken over T. Varying h{9), it 
follows that for all 9 we have 



o^'^qaixi, r, 9) dxi 


dr 


0 


and similarly for qp. Since (r, 9) are polar normal coordinates in Mo, the curves r ha (r, 9) 
are geodesics in Mo. Denoting the expression in brackets by fa{£,9) and varying the 
center of polar normal coordinates in Theorem 13 .1 1 and varying 9, we obtain that 


/„( 7 (r))exp 


- Xds 


dr = 0 


for all geodesics 7 in Mo which begin and end at points of dMo- This shows the vanishing 
of the geodesic ray transform of the function jd with constant attenuation —A. Since this 
transform is injective on simple two-dimensional manifolds |SU11] . we have /„ = 0 for all 
A. Thus 

/*oo 

/ e^^^'^qaixi, r, 9) dxi = 0 


for all A > 0, r and 9. Uniqueness for the Fourier transform in xi shows that go = 0 in 
M. We obtain g^ = 0 in M by the exact same argument. 

We have arrived at the following two equations in M: 

-lA(ai - 0:2) - \{d{ai + 02), d(ai - ^2)) + - £2/12) = 0 , 

— 1 A(/ 3 i — 132) — \{d{l 3 i + 132), d(/ 3 i — 132)) + — £2/12) = 0 . 
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Let u = {si/e 2 Y^‘^ and v = Then |(ai — 0 ^ 2 ) = log u and \{Yi — (32) = log n, 

and the eqnations become 

-A(log u) - (£l£ 2 )"^^^(d(£l£ 2 )^^^ d(log u)) + - £ 2 /^ 2 ) = 0, 

- A(log v) - d(log v)) + - £ 2 /^ 2 ) = 0. 

Multiplying the hrst equation by (£ 162 )^'^^ and the second by (/ii/i 2 )^'^^ and using that 
6{adw) = —aAw — {da,dw), we obtain the equations 

5{{eie2Y^‘^d{\ogu)) + u^{eie2Y^'^£2l^2{u‘^v^ - 1 ) = 0 , 

(5((/ii/i2)^'^^d(logn)) +a;^(/ii/i2)^/^£2h2(MV - 1 ) = 0. 

Since d{\ogu) = u~^du, d(logn) = v~^dv, we see that u and v satisfy the semilinear 
elliptic system 

S{e2du) + ijj'^e\iJi2{u^v'^ — 1 )m = 0 , 

6{fi2dv) + a;^£ 2 /^ 2 ('^^'^^ “ 1 )"^ = 0 - 

The assumptions on Sj and fij on dM imply that u = v = 1 and di^u = d^v = 0 on dM. 
Also, the above equations imply that the pair {u, v) = ( 1 , 1 ) is a solution of the semilinear 
system in all of M. Unique continuation holds for this system (see for instance |KSUn9( 
Appendix B]), and we obtain u = 1 and n = 1 in M. This proves that £1 = £2 and 
/ii = /i 2 in M as required. □ 

Proof of Theorem M.^ Assume for simplicity that Q C {xs > 0} and that Xq = 0. Let 
ip{x) = log |x| be an LCW in As explained in |KS13( Section 3B], if one chooses new 
coordinates 

X 

yi = log |x|, y' = -ri 
\x\ 

where i/' G = {x G ; X 3 > 0}, then the LCW becomes ip{y) = yi and the Euclidean 
metric becomes 

i 90(9') ) 

where go is the standard metric on the sphere Since manifolds {x G ; X 3 > c > 0} 
are simple, we have reduced matters to Theorem ll.il This proves Theorem 11.21 □ 

5. Appendix 

For a compact oriented smooth manifold M with smooth boundary, define 
H^ = {ue L‘^{M, AM) I Am G L^M, AM)}, 
where A is the Hodge Laplacian. We let 

In this appendix, we will prove the following trace theorem for 

Lemma 5.1. The relative and absolute trace maps u 1 —)■ {tu,tdu) and m 1 —)■ {t * u,t6 * u) 
initially defined on the space AM) have extensions as hounded linear maps from 

LTa to H-^^^{dM,AdM) x H-^/^{dM,AdM). The set A°°(M, AM) is dense in and 
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thus these extensions are unique. The trace maps u i—)■ (tw, tdu) and u ^ (t * u,td * u) 
have extensions to the same spaces, and one has for any u G 

(5.1) t{du) = dQuitu), 

(5.2) t{d * u) = ddM(t * u). 

Lemma [2.21 follows immediately from this result. To begin the proof of Lemma [5.11 we 
will need the following two lemmas. Here u is meant to signify the normal vector held on 
dM, extended into M in a standard way. 

Lemma 5.2. Let j > 0. For any v G {M, AM) and f,h ^ H^^^/'^{dM,AdM), there 
is a unique u G AM) that solves 

Am = V in M, tu = f, t * u = h. 

One has ||M||j^j +2 < C{\\v\\Hj + WfWw+s/^ + WhUnj+a/^)■ 

Proof. This is |Sc951 Theorem 3.4.10]. □ 

Lemma 5.3. For any {f,h) E H^{dM,AdM)xH^{dM,AdM), there is v E {M, AM) 
such that 

1^11^2 < ||/||rf3/2 + ||h||^l/2 

and tiyV = f, ti^, * v = tiyd * n = 0 , and ti^dv = h. 


Proof. Let hrst / = 0. Any boundary point has a neighborhood U (Z M where boundary 
normal coordinates (xi,... ,x„) exist, so that the boundary is given by = 0} and x„ 
coincides with the direction of u. We hrst dehne v in U using these coordinates. Then v 
is represented as a sum of objects of the form 


Vn,...,i^dxi,^ A ... A dx4 


for 1 < ifc < n, where are functions. Moreover, h can be represented as a sum of 

objects of the form 

A ... A dxi,^, 

for I < ik < n. 

We will dehne to be identically zero if any ik = n, and choose the remaining 

to be functions such that 




0 and 


on U r\ dM. Observe that when chosen this way all components of the form A 

... Adxjj, vanish on the boundary and iuVii,...,i^,dxi.^ A... Adx^ vanish in a neighborhood of 
the boundary. We can arrange for ^ \\hii,...,ik\\mZ{undM), as an inequality 

for functions. Using a suitable partition of unity we can dehne v near dM and extend 
smoothly as an element v E P[‘^{M, AM), and so we will obtain the inequality 

\\'v\\h^(M,AM) ^ \\h\\H^/^{aM,AdM) 

for V and h as forms. Note that our dehnition is made to guarantee that v\dM = 0, 
tV,yV = h, and iyV = 0 in a neighbhorhood of the boundary of M. Now by Lemma 3.4 
of [UbT13| . 

tiydv = tVvV\\ + Stv\\ — d'tiyV 
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where S' is a bounded endomorphism on AM and d! denotes the exterior derivative on 
dM. The last two terms vanish because v is zero on dM. 

Then we are left with 

tiydv = = tVuV = h. 

Moreover, ti^v = tiy * v = Q since v is zero on the boundary of M, and ti^d * u = 0 since 

tiyd *v = — d'tiy * v. 

and (*u)|| = *v±_ = 0 in a neighbourhood of the boundary. This proves the lemma when 
/ = 0. The next step is to consider the case where / is nonzero but h = 0. An argument 
in boundary normal coordinates similar to the one above yields a form v with ti^v = /, 
tiy *v = tiyd * V = tiudv = 0, and ||u||_f /2 < ||/||//3/2. Combining these two cases proves 
the lemma. □ 

Proof of Lemma \5.1[ Let u, u G AM). Then using the integration by parts formu¬ 

las fl2.2p and fl2.3p we get the identity 

(5.3) {u\Av)m — (Am|u)m = (u\iudv — v A 6v)qm + {ly A 6u — i^du\v)dM- 

We can rewrite the boundary terms in terms of the relative and absolute boundary values 
of M, by using the following formulas valid for any 1-form ^ and any fc-form rj (see for 
example |CST13j ): 

i^T) = (—i A *7], * * 7] = (— St] = (—* d*ri, 

{*u, *v) = {u, v), A u, v) = {u, i^v), {Au, v)\oM = {tiuU, tv) 

After a computation (hrst assuming that u and v are fc-forms and then summing over k), 
fl5.3p can be written as 

(5.4) {u\Av)m — (Am|u)m = {tu\ti^dv)aM + (t* u\ti^d * v)3m 

+ {t5u\tiyv)dM + {t5 * u\ti^ * v)oM- 

Now if M G Ha, we wish to define the relative boundary values of u by fl5.4p . Given 
any (/, h) G H^{dM, AdM) x H^{dM, AdM), we use LemmaESlto choose v = such 
that = /, = h, til, * V = ti„d * u = 0 and ||u||_h -2 ^ ||/||j:^3/2 -f ||h||j|^i/ 2 . By fl5.4p 

we have for any u G AM) 

itu\h)dM = {u\Avo^h)M - {Au\vo,h)M, 

{tSu\f)aM = (m|Au/,o)m - (Am|u/,o)m- 

Since 

||fM||j^-i /2 = sup \{tu\h)9M\, \\t6u\\H-3/2 = sup \{t6u\f)9M\, 

11^11^1/2=1 Il/llj:j3/2=l 

it follows that 

||fM||j:^-l/2 ||tdM||^-3/2 < ||m||h^, U G H^{M,AM). 

Thus for any u G Ha, we may use the above formulas to dehne {tu,t6u) as an element 
of H~^^‘^{dM, A{dM)) x H~^/‘^{dM, A{dM)). Applying this argument to *u shows that 
{t *u,t6* u) are well dehned for u G Ha in a similar way. The identity fl5.4p then remains 
true for any u G Ha and any w G H^{M, AM). 

We now claim that C^{M, AM) is dense in Ha- Given this fact, it follows that the 
maps u i-A {tu, t5u) and u ha {t*u,t6 *u) are bounded from Ha to x and are 
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uniquely defined by the corresponding maps acting on C°°. Finally, (15.Ih and (15.2^ are 
true for u G AM) by the fact that d commutes with the pull-back to the boundary 

t and thus they extend to u G -ffA by density. 

It remains to show that AM) is dense in Ha- Suppose that u G Ha- Then 

Am G L2(M,AM) and tu,t*u e H-^/‘^{dM, AdM), so for any e > 0 there exist v G Vt{M) 
and f,hG Vt{dM) such that 

(5.5) II Am — M||i 2 -I- 11/ — tu\\H-i /2 + \\h — t * u\\jj-i /2 < e- 

Now by Lemma 15.21 there exists a smooth form u' such that Am' = v, with tu' = f 
and t * u' = h. Then if w G H‘^{M, AM) satishes tw = t * w = 0, which implies 
tiyW = * w = Q, we get by fl5.4p that 

(m — u'\/Aw)m = (A(m — m ')| m;)m + {t{u — u')\tii,dw)dM + — u')\tii,d * w)om- 

Combining this with inequality fl5.5p we obtain 

|(m — m'|Am;)m| < Ce\\w\\H^- 

Now by Lemma 15.21 every form z G L‘^{M, AM) can be written as Aw for some w with 
tw = t *w = 0- Moreover, we will have 


w\\h^ ^ Iklk^. 


Thus we get that 


M — M \\l2 


< 

r\j 


e 


and hence 

\\u-u'\\h^ <£- 

Therefore there is a sequence of smooth forms Uk such that m^ —)■ m in Ha- 


□ 
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